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Abstract 

Sufficient conditions are developed, under which the compound Poisson distribution has maximal en- 
tropy within a natural class of probability measures on the nonnegative integers. Recently, one of the 
authors [O. Johnson, Stoch. Proc. AppL, 2007] used a semigroup approach to show that the Poisson has 
maximal entropy among all ultra-log-concave distributions with fixed mean. We show via a non-trivial 
extension of this semigroup approach that the natural analog of the Poisson maximum entropy property 
remains valid if the compound Poisson distributions under consideration are log-concave, but that it 
fails in general. A parallel maximum entropy result is established for the family of compound binomial 
measures. Sufficient conditions for compound distributions to be log-concave are discussed and appli- 
cations to combinatorics are examined; new bounds are derived on the entropy of the cardinality of 
a random independent set in a claw-free graph, and a connection is drawn to Mason's conjecture for 
matroids. The present results are primarily motivated by the desire to provide an information-theoretic 
foundation for compound Poisson approximation and associated limit theorems, analogous to the cor- 
responding developments for the central limit theorem and for Poisson approximation. Our results 
also demonstrate new links between some probabilistic methods and the combinatorial notions of log- 
concavity and ultra-log-concavity, and they add to the growing body of work exploring the applications 
of maximum entropy characterizations to problems in discrete mathematics. 
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1. Introduction 



The primary motivation for this work is the development of an information-theoretic approach to discrete 
limit laws, specifically those corresponding to compound Poisson limits. Recall that the classical central 
limit theorem can be phrased as follows: If X\, X 2 , ■ ■ ■ are independent and identically distributed, 
continuous random variables with zero mean and unit variance, then the entropy of their normalized 
partial sums S n = ^= increases with n to the entropy of the standard normal distribution, 

which is maximal among all random variables with zero mean and unit variance. More precisely, if f n 
denotes the density of S n and (f> the standard normal density, then, as n — > oo, 

h>(fn) t h{4>) = sup{h(f ) : densities / with mean and variance 1}, (1) 

where h(f) = — J f log / denotes the differential entropy, and 'log' denotes the natural logarithm. 
Precise conditions under which (JXJ) holds are given by Barron [B[ and Artstein et al. [2]; also see 



27, 56, 44] and the references therein. 



Part of the appeal of this formalization of the central limit theorem comes from its analogy to the 
second law of thermodynamics: The "state" (meaning the distribution) of the random variables S n 
evolves monotonically, until the maximum entropy state, the standard normal distribution, is reached. 
Moreover, the introduction of information-theoretic ideas and techniques in connection with the entropy 
has motivated numerous related results (and their proofs), generalizing and strengthening the central 
limit theorem in different directions; see the references above for details. 

Recently, some discrete limit laws have been examined in a similar vein, but, as the discrete entropy 
H(P) = —J2xP( x )\°g-P( x ) f° r probability mass functions P on a countable set naturally replaces 
the differential entropy h(f ), many of the relevant analytical tools become unavailable. For Poisson 
convergence theorems, of which the binomial-to-Poisson is the prototypical example, an analogous 
program has been carried out in 5^, 23, 41, 28[. Like with the central limit theorem, there are two 



aspects to this theory - the Poisson distribution is first identified as that which has maximal entropy 
within a natural class of probability measures, and then convergence of appropriate sums of random 
variables to the Poisson is established in the sense of relative entropy (or better still, approximation 
bounds are obtained that quantify the rate of this convergence) . 

One of the main goals of this work is to establish a starting point for developing an information-theoretic 
framework for the much more general class of compound Poisson limit theoremsil To that end, our first 
main result, given in Section 2, provides a maximum entropy characterization of compound Poisson 
laws, generalizing Johnson's characterization [28| of the Poisson distribution. It states that if one looks 
at the class of all ultra-log-concave distributions on Z + with a fixed mean, and then compounds each 
distribution in this class using a given probability measure on N, then the compound Poisson has 
maximal entropy in this class, provided it is log-concave. 



3 Recall that the compound Poisson distributions are the only infinitely divisible distributions on , and also they are 
(discrete) stable laws [5511 . In the way of motivation, we may also recall the remark of Gnedenko and Korolev [Tsl pp. 211- 
215] that "there should be mathematical . . . probabilistic models of the universal principle of non-decrease of uncertainty," 
and their proposal that we should "find conditions under which certain limit laws appearing in limit theorems of probability 
theory possess extremal entropy properties. Immediate candidates to be subjected to such analysis are, of course, stable 
laws." 
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Having established conditions under which a compound Poisson measure has maximum entropy, in a 
companion work Q we consider the problem of establishing compound Poisson limit theorems as well 
as finite-n approximation bounds, in relative entropy and total variation. The tools developed in the 
present work, and in particular the definition and analysis of a new score function in Section 3, play a 
crucial role in these compound Poisson approximation results. 

In a different direction, in Section 6 we demonstrate how the present results provide new links between 
classical probabilistic methods and the combinatorial notions of log-concavity and ultra-log-concavit; 
Log-concave seq uences are well-studied objects in combinatorics, see, e.g., the surveys by Brenti 
and Stanley [54| . Additional motivation in recent years has come from the search for a theory of 
negative dependence. Specifically, as argued by Pemantle a theory of negative dependence has long 
been desired in probability and statistical physics, in analogy with the theory of positive dependence 
exemplified by the Fortuin-Kasteleyn-Ginibre (FKG) inequality [l8| (earlier versions were developed 
by Harris [25| and, in combinatorics, by Kleitman [38]). But the development of such a theory is 
believed to be difficult and delicate. For instance, Wagner (HtJ recently formalized a folklore conjecture 
in probability theory (called by him the "Big Conjecture") which asserts that, if a probability measure 
on the Boolean hypercube satisfies certain negative correlation conditions, then the sequence {pt} of 
probabilities of picking a set of size k, is ultra-log-concave. This is closely related to Mason's conjecture 
for independent sets in matroids, which asserts that the sequence {Ik}, where Ik is the number of 
independent sets of size k in a matroid on a finite ground set, is ultra-log-concave. Soon after, the 
"Big Conjecture" was falsified, by Borcea, Branden and Liggett 0] and by Kahn and Neiman [3^] , who 
independently produced counterexamples. In the other direction, very recently (while we were revising 
this paper), Lenz [42] proved a weak version of Mason's conjecture. In Section 6 we describe some 
simple consequences of our results in the context of matroid theory, and we also discuss an application 
to bounding the entropy of the size of a random independent set in a claw-free graph. 

Before stating our main results in detail, we briefly mention how this line of work connects with the 
growing body of work exploring applications of maximum entropy characterizations to discrete math- 
ematics. The simplest maximum entropy result states that, among all probability distributions on a 
finite set 5, the uniform has maximal entropy, log 1 51. While mathematically trivial, this result, com- 
bined with appropriate structural assumptions and various entropy inequalities, has been employed as 
a powerful tool and has seen varied applications in combinatorics. Examples include Radhakrishnan's 



entropy-based proof [5l| of Bregman's theorem on the permanent of a 0-1 matrix, Kahn's proof [3l| of 
the result of Kleitman and Markowski (35| on Dedekind's problem concerning the number of antichains 
in the Boolean lattice, the study by Brightwell and Tetali Q of the number of linear extensions of 
the Boolean lattice (partially confirming a conjecture of Sha and Kleitman [13]), and the resolution 
of several conjectures of Imre Ruzsa in additive combinatorics by Madiman, Marcus and Tetali [45^ . 
However, so far, a limitation of this line of work has been that it can only handle problems on finite sets. 
As modern combinatorics explores more and more properties of countable structures - such as infinite 
graphs or posets - it is natural that analogues of useful tools such as maximum entropy characterizations 
in countable settings should develop in parallel. It is particularly natural to develop these in connection 
with the Poisson and compound Poisson laws, which arise naturally in probabilistic combinatorics; see, 
e.g., Penrose's work [H3] on geometric random graphs. 

Section 2 contains our two main results: The maximum entropy characterization of log-concave com- 
pound Poisson distributions, and an analogous result for compound binomials. Sections 3 and 4, 
respectively, provide their proofs. Section 5 discusses conditions for log-concavity, and gives some ad- 
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ditional results. Section 6 discusses applications to classical combinatorics, including graph theory and 
matroid theory. Section 7 contains some concluding remarks, a brief description of potential general- 
ization and extensions, and a discussion of recent, subsequent work by Y. Yu [59], which was motivated 
by preliminary versions of some of the present results. 



2. Maximum Entropy Results 

First we review the maximum entropy property of the Poisson distribution. 

Definition 2.1. For any parameter vector p = (pi,p2, . . . ,p n ) with each pi G [0,1], the sum of inde- 
pendent Bernoulli random variables Bi ~ Bern(pi), 

n 

S n = ^ Bi, 

i=l 

is called a Bernoulli sum, and its probability mass function is denoted by b p (x) := Pr{S n ~ x}, for 
x = 0, 1, . . .. Further, for each A > 0, we define the following sets of parameter vectors: 

Pn(A) = {pe[0,l] n : Pi+P2 + ---+p n =X} and P^X) = (J P n (X). 

n>l 



Shepp and Olkin 53] showed that, for fixed n > 1, the Bernoulli sum 6 P which has maximal entropy 



among all Bernoulli sums with mean A, is Bin(n, X/n), the binomial with parameters n and X/n, 

Jf(Bin(n,A/n)) =max{ff(6 p ) : peP„(A)}, (2) 

where H(P) = — P( x ) 1°S -P( x ) denotes the discrete entropy function. Noting that the binomial 
Bin(n, X/n) converges to the Poisson distribution Po(A) as n — > oo, and that the classes of Bernoulli 
sums in are nested, {6 P : p G P„(A)} G {b p : p G V n+ i(X)}, Harremoes noticed that a simple 
limiting argument gives the following maximum entropy property for the Poisson distribution: 

H(Po(X)) = sup (6 p ) : p G 7>oo(A)}. (3) 

A key property in generalizing and understanding this maximum entropy property further is that 
of ultra-log-concavity; cf. [49] . The distribution P of a random variable X is ultra-log-concave if 
P(x)/H\(x) is log-concave, that is, if, 

xP(x) 2 > (x + l)P(x + l)P(x - 1), for all x > 1. (4) 

Note that the Poisson distribution as well as all Bernoulli sums are ultra-log-concave. A non-trivial 
property of the class of ultra-log-concave distributions, conjectured by Pemantle and proved by 



Liggett [43J (cf. Gurvits [20| and Kahn and Neiman (33J), is that it is closed under convolution. 

Johnson [2q | recently proved the following maximum entropy property for the Poisson distribution, 
generalizing ([3]): 

iJ(Po(A)) = max | H{P) : ultra- log-concave P with mean a|. (5) 
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As discussed in the Introduction, we wish to generalize the maximum entropy properties © and (O to 
the case of compound Poisson distributions on Z + . We begin with some definitions: 

Definition 2.2. Let P be an arbitrary distribution on Z + = {0, 1, . . .}, and Q a distribution on N = 
{I, 2, . . .}. The Q-compound distribution CqP is the distribution of the random sum, 

XX ( 6 ) 

where Y has distribution P and the random variables {Xj} are independent and identically distributed 
(i.i.d.) with common distribution Q and independent ofY. The distribution Q is called a compounding 
distribution, and the map P n- CqP is the Q-compounding operation. The Q-compound distribution 
CqP can be explicitly written as the mixture, 

oo 

c Q p(x) = Y. p (y)Q* y ( x )> x ^°> ( 7 ) 

y=o 

where Q* 3 (x) is the jth convolution power of Q and Q*° is the point mass at x = 0. 

Above and throughout the paper, the empty sum Ylj=i(' ' ' ) i s taken to be zero; all random variables 
considered are supported on Z + = {0, 1, . . .}; and all compounding distributions Q are supported on 

n \\. ■>....}. 

Example 2.3. Let Q be an arbitrary distribution on N. 

1. For any < p < 1, the compound Bernoulli distribution CBern (p, Q) is the distribution of the 
product BX , where B ~ Bern(p) and X ~ Q are independent. Lt has probability mass function 
CqP, where P is the Bern (p) mass function, so that, CqP(0) = 1 — p and CqP(x) = pQ(x) for 
x>l. 

2. A compound Bernoulli sum is a sum of independent compound Bernoulli random variables, all 
with respect to the same compounding distribution Q: Let X\, X2, ■ ■ ■ , X n be i.i.d. with common 
distribution Q and B\, Bi, ■ ■ ■ , B n be independent Bern(pi). We call, 

/jBjXj = 22 Xj, 
i=\ j=l 

a compound Bernoulli sum; in view of its distribution is CQb p , where p = (pi,p2, ■ ■ • ,Pn)- 

3. Ln the special case of a compound Bernoulli sum with all its parameters pi — p for a fixed p € [0, 1], 
we say that it has a compound binomial distribution, denoted by CBin(n,p, Q). 

4. Let Tl\(x) = e- x X x /x\, x > 0, denote the Po(A) mass function. Then, for any A > 0, the 
compound Poisson distribution CPo(A, Q) is the distribution with mass function CqTI\: 

00 00 e~ x \i 

c Q n A (x) = ^n A (^(x)^— ^w, x>o. (8) 
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In view of the Shepp-Olkin maximum entropy property @ for the binomial distribution, a first natural 
conjecture might be that the compound binomial has maximum entropy among all compound Bernoulli 
sums Cq&p with a fixed mean; that is, 

£T(CBin(n,A/n,Q)) = max {iJ(C Q o p ) : P eP„(A)}. (9) 

But, perhaps somewhat surprisingly, as Zhiyi Chi has noted, © fails in general. For example, 
taking Q to be the uniform distribution on {1, 2}, p = (0.00125, 0.00875) and A — p\ + P2 = 0.01, direct 
computation shows that, 

F(CBin(2, A/2, Q)) < 0.090798 < 0.090804 < H(C Q b p ). (10) 

As the Shcpp-Olkin result (HJ) was only seen as an intermediate step in proving the maximum entropy 
property of the Poisson distribution §5$ , we may still hope that the corresponding result remains true 
for compound Poisson measures, namely that, 

Jf(CPo(A,Q)) = sup{H(C Q b p ) : p e Poo(A)}. (11) 

Again, (jlll) fails in general. For example, taking the same Q, A and p as above, yields, 

ff(CPo(A,Q)) < 0.090765 < 0.090804 < H(C Q b p ). 

The main purpose of the present work is to show that, despite these negative results, it is possible 
to provide natural, broad sufficient conditions, under which the compound binomial and compound 
Poisson distributions can be shown to have maximal entropy in an appropriate class of measures. 

Our first result (a more general version of which is proved in Section 3) states that, as long as Q and 
the compound Poisson measure CPo(A, Q) are log-concave, the maximum entropy statement analogous 
to ([3]) remains valid in the compound Poisson case: 

Theorem 2.4. If the distribution Q on N and the compound Poisson distribution CPo(A, Q) are both 
log- concave, then, 

iJ(CPo(A, Q)) = max ^H{CqP) : ultra-log- concave P with mean a|. 

The notion of log-concavity is central in the development of the ideas in this work. Recall that the 
distribution P of a random variable X on Z + is log-concave if its support is a (possibly infinite) interval 
of successive integers in Z+, and, 

P{x) 2 > P(x + l)P(x- 1), for all x > 1. (12) 

We also recall that most of the commonly used distributions appearing in applications (e.g., the Pois- 
son, binomial, geometric, negative binomial, hypergeometric logarithmic series, or Polya-Eggenberger 
distribution) are log-concave. 

Note that ultra-log-concavity of P, defined as in equation is more restrictive than log-concavity, 
and it is equivalent to the requirement that ratio P/H\ is a log-concave sequence for some (hence all) 
A > 0. 

Our second result states that (j9]) does hold, under certain conditions on Q and CBin(n, A, Q): 
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Theorem 2.5. If the distribution Q on N and the compound binomial distribution CBin(n, X/n, Q) are 
both log-concave, then, 

H{CBm{n,X/n,Q)) =max( H(C Q b p ) : p e V n (X)\, 

as long as the tail of Q satisfies either one of the following properties: (a) Q has finite support; or (b) Q 
has tails heavy enough so that, for some p, j3 > and N > 1, we have, Q(x) > p x , for all x > N . 



The proof of Theorem 12.51 is given in Section 01 As can be seen there, conditions (a) and (b) are 
introduced purely for technical reasons, and can probably be significantly relaxed; see Section 7 for a 
further discussion. 

It remains an open question to give necessary and sufficient conditions on A and Q for the compound 
Poisson and compound binomial distributions to have maximal entropy within an appropriately defined 
class. As a first step, one may ask for natural conditions that imply that a compound binomial or 
compound Poisson distribution is log-concave. We discuss several such conditions in Section 5. 

In particular, the discussion in Section 5 implies the following explicit maximum entropy statements. 

Example 2.6. 1. Let Q be an arbitrary log-concave distribution on N. Then Lemma \5. 1\ combined 
with Theorem \2.5\ implies that the maximum entropy property of the compound binomial distribu- 
tion in equation (|9]) holds, for all X large enough. That is, the compound binomial CBin(n, X/n, Q) 
has maximal entropy among all compound Bernoulli sums Cq6 p with p\ + p2 + • • • + p n = X, as 
long asX> Q(1 "?j^ (2) ■ 

2. LetQ be an arbitrary log-concave distribution onN. Then Theorem ] 5. 5\ combined with Theorem \2.4\ 
implies that the maximum entropy property of the compound Poisson CPo{X, Q) holds if and only 
if A > 

V A - Q(l) 2 ■ 

As mentioned in the introduction, the above results can be used in order to gain better understanding 
of ultra-log-concave sequences in combinatorics. Specifically, as discussed in more detail in Section 6, 
they can be used to estimate how "spread out" these sequences in terms of the entropy. 



3. Maximum Entropy Property of the Compound Poisson Distribution 



Here we show that, if Q and the compound Poisson distribution CPo(A, Q) = CqI1\ are both log- 
concave, then CPo(A, Q) has maximum entropy among all distributions of the form CqP, when P has 
mean A and is ultra- log-concave. Our approach is an extension of the 'semigroup' arguments of [28j . 

We begin by recording some basic properties of log-concave and ultra-log-concave distributions: 

(i) If P is ultra-log-concave, then from the definitions it is immediate that P is log-concave. 

(ii) If Q is log-concave, then it has finite moments of all orders; see [3(3, Theorem 7]. 

(m) If X is a random variable with ultra- log-concave distribution P, then (by (i) and (ii)) it has finite 
moments of all orders. Moreover, considering the covariance between the decreasing function 
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P(x + l)(x + 1)/P(x) and the increasing function x(x — 1) • • • (x — n), shows that the falling 
factorial moments of P satisfy, 

E[{X) n ] := E[X{X - 1) • • ■ (X - n + 1)] < (E(X)) n ; 



see [28| and [24| for details. 
(iv) The Poisson distribution and all Bernoulli sums are ultra-log-concave. 

Recall the following definition from [28| : 

Definition 3.1. Given a € [0,1] and a random variable X ~ P on Z + with mean A > 0, let U a P 
denote the distribution of the random variable, 

x 

^ Bi + Zx(l- a ), 
i=l 

where the Bi are i.i.d. Bern (a), ^A(i-a) has distribution Po(A(l — a)), and all random variables are 
independent of each other and of X . 

Note that, if X ~ P has mean A, then U a P has the same mean. Also, recall the following useful relation 
that was established in Proposition 3.6 of [28[: For all y > 0, 

^-U a P(y) = - (X(U a P(y) - U a P(y - 1) - ((y + l)U a P(y + 1) - yU a P(y))) . (13) 



Next we define another transformation of probability distributions P on Z 



+ ■ 



Definition 3.2. Given a G [0, 1], a distribution P on 1i + and a compounding distribution Q on N, let 
U§P denote the distribution CQU a P: 

oo 

UgP(x) := C Q U a P(x) = ]T U a P(y)Q*y(x), x>0. 

v=o 



Work of Chafa'i [Tl| suggests that the semigroup of Definition 13.11 may be viewed as the action of the 
M/M /oo queue. Similarly the semigroup of Definition 13.21 corresponds to the marginal distributions of 
a continuous-time hidden Markov process, where the underlying Markov process is the M/M /oo queue 
and the output at each time is obtained by a compounding operation. 

Definition 3.3. For a distribution P on Z + with mean v, its size-biased distribution P# on Z + is 
defined by 

p#{y) = (y + my + V 

V 

An important observation that will be at the heart of the proof of Theorem 12 .41 below is that, for a = 0, 
UqP is simply the compound Poisson measure CP (A, Q), while for a = 1, P = CqP. The following 
lemma gives a rough bound on the third moment of U®P: 
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Lemma 3.4. Suppose P is an ultra-log-concave distribution with mean A > on %+, and let Q be a 
log-concave compounding distribution on N. For each a € [0,1], let W a ,V a be random variables with 
distributions U®P = CqIIcP and CQ(U a P)# , respectively. Then the third moments E(W 3 ) and E(V 3 ) 
are both bounded above by, 

\q 3 + 3\ 2 qi q 2 + \ 3 ql 
where 171,(72,93 denote the first, second and third moments of Q, respectively. 

Proof. Recall that, as stated in properties (ii) and (Hi) in the beginning of Section [3J Q has finite 
moments of all orders, and that the nth falling factorial moment of any ultra-log-concave random 
variable Y with distribution R on Z + is bounded above by (E(Y)) n . Now for an arbitrary ultra-log- 
concave distribution R, define random variables Y ~ R and Z ~ CqR. If ri, r 2 , r% denote the first three 
moments of Y ~ R } then, 

E{Z 3 ) = q 3 r 1 +3q 1 q 2 E[(Y) 2 }+qlE[(Y) 3 } 

< q3ri + 3qiq 2 rt + qM. (14) 

Since the map U a preserves ultra-log-concavity [28}, if P is ultra-log-concave then so is R — U a P, so 
that (fT4")) gives the required bound for the third moment of W a , upon noting that the mean of the 
distribution U a P is equal to A. 

Similarly, size-biasing preserves ultra- log-concavity; that is, if R is ultra- log-concave, then so is R#, 
since R*{x + l){x + l)/R#{x) = (R(x+ 2)(x+ 2)(x+ l))/(R(x + l)(x + 1)) = R(x + 2)(x + 2)/R{x + 1) 
is also decreasing. Hence, R' = (U a P)# is ultra-log-concave, and (fT4")l applies in this case as well. In 
particular, noting that the mean of Y' ~ R' = (U a P)# = R# can be bounded in terms of the mean of 
Y - R as, 

E(Y') - y, (£+l)gag(£+l) = E K Y h] < >L = A 
* A E ( Y ^ ~ A 

the bound (|14p yields the required bound for the third moment of V a . □ 

In p8| . the characterization of the Poisson as a maximum entropy distribution was proved through 
the decrease of its score function. In an analogous way, we define the score function of a Q-compound 
random variable as follows, cf. 

Definition 3.5. Given a distribution P on Z + with mean X, the corresponding Q-compound distribution 
CqP has score function defined by, 

C Q (P*)(x) 
ri ^ pW= C Q P(x) - 1 ' X ^°- 



More explicitly, one can write 



ri - CQPix) = "x^moW) 1 (15) 

Notice that the mean of of r\^ Q p with respect to CqP is zero, and that if P ~ Po(A) then r\^c Q p{x) = 0. 
Further, when Q is the point mass at 1 this score function reduces to the "scaled score function" 
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introduced in 41 1. But, unlike the scaled score function and an alternative score function given in QJ, 
this score function is not only a function of the compound distribution CqP, but also explicitly depends 
on P. A projection identity and other properties of r\ t c Q P are proved in [4|. 

Next we show that, if Q is log-concave and P is ultra-log-concave, then the score function rx,c Q p{%) is 
decreasing in x. 

Lemma 3.6. If P is ultra-log- concave and the compounding distribution Q is log-concave, then the 
score function ri^c Q p{x) of CqP is decreasing in x. 

Proof. First we recall Theorem 2.1 of Keilson and Sumita (37[, which implies that, if Q is log-concave, 
then for any m > n, and for any x: 

Q* m (x + l)Q* n {x) - Q* m (x)Q* n {x + 1) > 0. (16) 

[This can be proved by considering Q* m as the convolution of Q* n and Q*( m - n ) j and writing 

Q* m (x + l)Q* n (x) - Q* m (x)Q* n (x + 1) 

= Q* (m ~ n) (l) (Q* n (x + 1 - l)Q* n {x) - Q* n (x - l)Q* n (x + 1) 



Since Q is log-concave, then so is Q* n , cf. [3J], so the ratio Q* n (x + 1)/Q* n (x) is decreasing in x, and 
(fTBT) follows.! 



By definition, ri,c Q p(%) > ^i,c q p{x + 1) if and only if, 

(j2(y + l)P{y + l)Q*y{x + l^j P^)Q* Z ^ 

= ^(y+ljPCy+W^IQ^W^ + lJ-Q^^ + lJQ**^)]. (17) 

y,z 

Noting that for y — z the term in square brackets in the double sum becomes zero, and swapping the 
values of y and z in the range y > z, the double sum in (fl7)) becomes, 

J2l(y + i)P{y + i)P(z) - (z + i)P(* + i)P(y)] [Q* y (x)Q* z (x + 1) - Q*y{ x + i)Q**(*)] . 

y<z 

By the ultra- log-concavity of P, the first square bracket is positive for y < z, and by equation (|16|) the 
second square bracket is also positive for y < z. □ 

We remark that, under the same assumptions, and using a very similar argument, an analogous result 
holds for some alternative score functions recently introduced in Q and in related work. 



Combining Lemmas 13.61 and 13.41 with equation (1131) we deduce the following result, which is the main 
technical step in the proof of Theorem 12.41 below. 
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Proposition 3.7. Let P be an ultra-log-concave distribution on Z + with mean A > 0, and assume that 
Q and CPo(A, Q) are both log-concave. Let W a be a random variable with distribution U§ P, and define, 
for all a € [0, 1], the function, 

E{a) i= E[-\QgC Q U x {W a )]. 

Then E(pt) is continuous for all a € [0, 1], it is differ entiable for a £ (0, 1), and, moreover, E'(a) < 
for a £ (0, 1). Ln particular, E(0) > E(l). 

Proof. Recall that, 

OO X 

UgP(x) = C Q U a P(x) = J2 U a P(y)Q*y(x) = ]T U a P(y)<?"(x), 

V=0 y=0 

where the last sum is restricted to the range < y < x, because Q is supported on N. Therefore, since 
U a P(x) is continuous in a [28], so is U^P(x), and to show that E{a) is continuous it suffices to show 
that the series, 

oo 

E(a) := E[-logC Q Il x (W a )} = - £ U%P{x) logC Q Il x (x), (18) 

x=0 

converges uniformly. To that end, first observe that log-concavity of CqIIa implies that Q(l) is nonzero. 
[Otherwise, if i > 1 be the smallest integer i such that Q(i) ^ 0, then CqIIa^ + I) = 0, but CQH\(i) and 
CqH\{2i) are both strictly positive, contradicting the log-concavity of CqIIa.] Since Q(l) is nonzero, 
we can bound the compound Poisson probabilities as, 

1 > C Q Ii x {x) = Y^[e~ X \ y /yl]Q* y (x) > e- x [X x /xl]Q(l) x , for all x > 1, 
v 

so that the summands in (|18[) can be bounded, 

< -logC Q n A (a;) < A + logx! -xlog(AQ(l)) < Cx 2 , x>l, (19) 

for a constant C > that depends only on A and Q(l). Therefore, for any N > 1, the tail of the series 
(fT8| can be bounded, 

oo „ 

< - uSP(x)logC Q Tl x (x) < CE[W 2 a l {Wa > N} ] < -E[Wl], 

x=N 

and, in view of Lemma 13.41 it converges uniformly. 

Therefore, E(a) is continuous in a, and, in particular, convergent for all a € [0, 1]. To prove that it is 
differentiable at each a € (0, 1) we need to establish that: (i) the summands in (|18l) are continuously 
differentiable in a for each x; and (ii) the series of derivatives converges uniformly. 

Since, as noted above, U®P(x) is defined by a finite sum, we can differentiate with respect to a under 
the sum, to obtain, 

lL U « P{X) = ~L CQUaP{x) = £ -^U a P(y)Q*y(x). (20) 

y=0 
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And since U a P is continuously differentiable in a E (0, 1) for each x (cf. [28|, Proposition 3.6] or equation 
(p~3j) above), so are the summands in (JTSj) , establishing (i); in fact, they are infinitely differentiable, which 
can be seen by repeated applications of (|13[) . To show that the series of derivatives converges uniformly, 
let a be restricted in an arbitrary open interval (e, 1) for some e > 0. The relation (|13p combined with 
(|2H| yields, for any x, 



d_ 

da 



USP(x) 



y=o 



E ( HU a P(y) - U a P(y - 1) - ((y + l)U a P(y + 1) - yU a P{y)) ) Q*^( 



= -- E ((» + 1 ) C7 « p (2/ + !) - AC7 «P(y)) (*) - Q^ +1 (x)) 
= -- E ((» + + 1) - AP«P(y)) 

+ E E ((» + l ) u <* p {y + x u a p(y)) q* v {* - «o 

Z / Q, Z / 

v=0 y=0 

« ^ A«?P(x) / 

a „ V \U£P(x - v) 



v=0 



A 
a 



(21) 



ti=0 



Also, for any x, by definition, 

|[/«P(x)r 1£/ Q p (x)| < C Q (E/ a P)#(x) + E/«P(x), 

where, for any distribution P, we write P^(y) = P(y + l)(y + 1)/A for its size-biased version. Hence 
for any N > 1, equations (l2~Tj) and (fH))) yield the bound, 

E T-^^WiogC Q n A ( a 



-N 



C\t 2 r x 1 

^ E ^{C'Q(^P) # (x) + ^P(x) + ^Q(i;)[C Q (C/ a P)#(x-i;) + C/«P(x- U )]} 

¥ a v=0 

-E\ Vl + Wl + X 2 + XV a + XW a )l { V a >N, W a >N, X>N} 



2C 



a 

a 



{E[V*l {Va > N} ] + E[WZl {Wa > N} ] + E[X 2 I {X > N} }} 
< j^{E[V*]+E\W*]+E[X a ]}, 



a 
C 
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where C, C > are appropriate finite constants, and the random variables V a ~ Cq(U P)# , W a ~ 
f/^P and X ^ Q are independent. Lemma 13.41 implies that this bound converges to zero uniformly 
in a £ (e, 1), as iV — >• oo. Since e > was arbitrary, this establishes that E(a) is differentiable for all 
a € (0, 1) and, in fact, that we can differentiate the series (|T5|> term-by-term, to obtain, 

= - J2^ U S P (^ogC Q U x (x) (22) 



n 5a 

a;— 



oo / a; \ 

2 UgP(x)r huSp (x) -Y,Q{v)U2P{x - v)r luSp {x - v) logC Q II A (x) 

oo / oo \ 

J2uSP(x)r liUSp (x) log<7 Q n A (x)-^Q(«)logCQn A (x + W ) , 



A °° 



x=0 \ v=0 



where the second equality follows from using (|2T|) above, and the rearrangement leading to the third 
equality follows by interchanging the order of (second) double summation and replacing x by x + v. 

Now we note that, exactly as in [28| . the last series above is the covariance between the (zero-mean) 
function r iu a p {x) and the function (log CqH\(x) — J2 V Q( v ) 1°S CqILx(x + y)), under the measure 
U§P. Since P is ultra-log-concave, so is U a P hence the score function r x u q p (x) is decreasing in 
x, by Lemma T3. 61 Also, the log-concavity of CqXV\ implies that the second function is increasing, and 
Chebyshev's rearrangement lemma implies that the covariance is less than or equal to zero, proving 
that E'(a) < 0, as claimed. 

Finally, the fact that E(Q) > E(l) is an immediate consequence of the continuity of E(a) on [0, 1] and 
the fact that E'(a) < for all a G (0, 1). □ 

Notice that, for the above proof to work, it is not necessary that CqIIa be log-concave; the weaker 
property that (logCgII\(x) — J2 V Qiy) log CqIIa(x + v)) be increasing is enough. 

We can now state and prove a slightly more general form of Theorem 12.41 Recall that the relative 
entropy between distributions P and Q on Z+, denoted by D(P\\Q), is defined by 

Theorem 3.8. Let P be an ultra-log- concave distribution on Z + with mean A. If the distribution Q on 
N and the compound Poisson distribution CqII\ are both log-concave, then 

D(c Q p\\c Q n x ) < H(c Q n x ) - h{c q p). 

Proof. As in Proposition [3Jl let W a - UgP = C Q U a P. Noting that W - CqU\ and Wi - CqP, we 
have 

H{C Q P) + D(C Q P\\C Q Il x ) = -SpogCon A (Wi)] 

< -E[logC Q n x (W Q )} 
= H(CqII\), 

where the inequality is exactly the statement that E(l) < E(0), proved in Proposition 13.71 □ 
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Since < D(C q P\\CqU x ), Theorem [2H immediately follows. 



4. Maximum Entropy Property of the Compound Binomial Distribution 

Here we prove the maximum entropy result for compound binomial random variables, Theorem l2.5l The 



proof, to some extent, parallels some of the arguments in 23] 47] 53], which rely on differentiating the 
compound-sum probabilities b p (x) for a given parameter vector p = (pi,p2, ■ ■ ■ ,Pn) (recall Definition ^. II 
in the Introduction) , with respect to an individual pi . Using the representation, 

a 
x=Q 

differentiating Cqb p (x) reduces to differentiating b p (x), and leads to an expression equivalent to that 
derived earlier in (12~T1) for the derivative of CnU a P with respect to a. 



Lemma 4.1. Given a parameter vector p = (j>i,P2, ■ ■ ■ ,p n ), with n > 2 and each < Pi < 1, let, 

f Pi + P2 . , Pi + P2 , 
Pt = I ^ ht ' 2 *>P3,---,Pn 

for t e [-(pi +pa)/2, (pi +Pa)/2]. I7ien, 



y=0 



71 

COM*) = (-2i)£&g(y) (Q^+ 2 )(x) -2Q^ +1 )(.x) + Q*^)) , (24) 



w/iere p = (p 3 , . . . ,p„). 

Proof. Note that the sum of the entries of pt is constant as t varies, and that pt = p for t = (pi —p2)/2, 
while p t = ({pi +p 2 )/2, (pi +p 2 )/2,_P3, • ■ • ,Pn) for t = 0. Writing k = pi +p 2 , b Pt can be expressed, 



'W,<» = [^-* a )Mw- 2 )+(*( 1 -|)+2/-)''p(//-ll 

2 



and its derivative with respect to i is 

9 



gt b pt (y) = -2t {b p (y - 2) - 2b p (y - 1) + 6 p (y)) . 
The expression (|23p for Cqb Pt shows that it is a finite linear combination of compound-sum probabilities 
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b Pt (x), so we can differentiate inside the sum to obtain, 

-^.c Q b Pt {x) = Qi b Pt(y)Q* v ( x ) 

v=o 



-2tJ2 {bp(y - 2) - 26p(y - 1) + b p (y)) Q*»(x) 
y=o 

n-2 

-2i^ 6g( V ) (g*^(x) - 2Q*^ +1 )(x) + Q*»(ao) , 



3/=0 

since 6p(y) = for y < — 1 and y > n — 1. □ 

Next we state and prove the equivalent of Proposition 13.71 above. Note that the distribution of a 
compound Bernoulli sum is invariant under permutations of the Bernoulli parameters pi. Therefore, 
the assumption p\ > pi is made below without loss of generality. 

Proposition 4.2. Suppose that the distribution Q on N and the compound binomial distribution 
CBin(n, X/n, Q) are both log-concave; let p = (pi,P2, • ■ • ,Pn) be a given parameter vector with n > 2, 
Pi + P2 + ■ • • + Pn = X > 0, and pi > pi; let Wt be a random variable with distribution CQb Pt ; and 
define, for all t £ [0, (pi — £>2)/2], the function, 

E(t) :=E[- logCq&p-ra, 

where p denotes the parameter vector with all entries equal to X/n. If Q satisfies either of the conditions: 
(a) Q finite support; or (b) Q has tails heavy enough so that, for some p, ft > and Nq > 1 , we have, 
Q(x) > p x , for all x > Nq, then E(t) is continuous for all t G [0, (p\ — f>2)/2], it is differentiable for t G 
(0, {pi -p 2 )/2), and, moreover, E'{t) < for t G (0, (p t -p 2 )/2). /n particular, E(0) > E(( Pl -p 2 )/2). 

Proof. The compound distribution Cqb Pt is defined by the finite sum, 

n 

and is, therefore, continuous in t. First, assume that Q has finite support. Then so does Cq6 p for any 
parameter vector p, and the continuity and differentiability of E(t) are trivial. In particular, the series 
defining E(t) is a finite sum, so we can differentiate term-by-term, to obtain, 

x— 
oo n— 2 

= 2^^&p( 2/ )(Q*fe+ 2 )( a; )-2g*fe+ 1 )( 2 :)+Q^(^))logCQM^) (25) 



n— 2 oo 

log C Q bp(z + v + w) - log C Q bp(z + v) 



y—0 z—0 v,w 



log C Q bp(z + w)+ log C Q bp(z) 



, (26) 
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where (f2"5| follows by Lemma [4.11 By assumption, the distribution Cgbp = CBin(n, A/n, Q) is log- 
concave, which implies that, for all z, v, w such that z + v + w is in the support of CBin(n, A/n, Q), 

Cq%0) < C Q bp(z + w) 



C Q bp(z + v) C Q bp(z + v + w)' 

Hence the term in square brackets in equation (1261) is negative, and the result follows. 

Now, suppose condition (b) holds on the tails of Q. First we note that the moments of Wt are all 
uniformly bounded in t: Indeed, for any 7 > 0, 

00 oo n n oo 

E\W?] =J2C Q b Pt (x)x'r = ^E^tfXr^ 7 < EX^O*)* 7 < C„g 7 , (27) 

x—0 x—0 y—0 y—0 x—0 

where C n is a constant depending only on n, and q 7 is the 7th moment of Q, which is of course finite; 
recall property (ii) in the beginning of Section [3] 

For the continuity of E(t), it suffices to show that the series, 

00 

E(t) := E[- log C Q bp(W t )} = -J2C Q b Pt (x)logC Q bv(x), (28) 

x=0 

converges uniformly. The tail assumption on Q implies that, for all x > No, 

n 

1 > Cgbpix) =Y,b?(y)Q* V (x) > A U - A/n) n - l Q{x) > A(l - A/n^V, 

y=0 

so that, 

0<-\ogC Q bp{x)<Cx p , (29) 



for an appropriate constant C > 0. Then, for N > No, the tail of the series (|28|) can be bounded, 
< - £ C Q b Pt (x) log Cqb^x) < CE[W t \ Wt > N} ] < ^E[wf +1 ] < ^C n q p+1 , 

x=N 



where the last inequality follows from (|27p . This obviously converges to zero, uniformly in t, therefore 
E(t) is continuous. 

For the differentiability of E(t), note that the summands in (|T5|) are continuously differentiable (by 
Lemma [4.ip . and that the series of derivatives converges uniformly in t; to see that, for N > No we 
apply Lemma T4. 1 1 together with the bound (j2l?|) to get, 

00 d 

Yl JT f C QKA x )^gC Q bp{3 



at 

x=N 



oo n 



< 2tJ2 I>p(2/) (Q* [v+2) {x) +2Q<y +1 \x) +Q* y (x)) Cx? 

x=N y=0 
n oo 

< 2C<X X [Q* {y+2 \x) +2Q<y +1 \x) +Q* y (xf) x?, 

y=0 X = N 
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which is again easily seen to converge to zero uniformly in t as N — > oo, since Q has finite moments of 
all orders. This establishes the differentiability of E(t) and justifies the term-by-term differentiation of 
the series (fT5)k the rest of the proof that E'(t) < is the same as in case (a). □ 

Note that, as with Proposition 13. 7\ the above proof only requires that the compound binomial distribu- 
tion CBin(n, X/n, Q) = CQbp satisfies a property weaker than log-concavity, namely that the function, 
logCQ&p(x) - J2v Q( v ) l°g CqM^ + v )> be increasing in x. 

Proof, (of Theorem I2.5[) Assume, without loss of generality, that n > 2. If p\ > p 2 , then Proposi- 
tion [O] says that, E((pi -p 2 )/2) < E(0), that is, 

OO OO 

-J2C Q b p (x)logC Q bp(x) < C Q b p .(x) log C Q bp(x), 

x=0 x=0 

where p* = ((|>i +p2)/2, (p± +P2)/2,P3, ■ ■ -Pn) and p = (X/n, . . . , X/n). Since the expression in the 
above right-hand-side is invariant under permutations of the elements of the parameter vectors, we 
deduce that it is maximized by pt = p. Therefore, using, as before, the nonnegativity of the relative 
entropy, 

H(C Q b p ) < ff(C Q & p ) + £(C Q 6 p ||C Q fep) 

OO 

= -^20 Q b p (x)logC Q bp(x) 

x=Q 
oo 

x=0 

- H(C Q bp) = H{CBm(n,X/n,Q)), 
as claimed. □ 

Clearly one can also state a slightly more general version of Theorem 12.51 analogous to Theorem 13.81 



5. Conditions for Log-Concavity 

Theorems 12.41 and 12.51 state that log-concavity is a sufficient condition for compound binomial and 
compound Poisson distributions to have maximal entropy within a natural class. In this section, we 
discuss when log-concavity holds. 

Recall that Steutel and van Harn [HBl Theorem 2.3] showed that, if {iQ(i)} is a decreasing sequence, then 
CPo(A, Q) is a unimodal distribution, which is a necessary condition for log-concavity. Interestingly, 
the same condition provides a dichotomy of results in compound Poisson approximation bounds as 
developed by Barbour, Chen and Loh [3(: If {iQ(i)} is decreasing, then the bounds are of the same 
form and order as in the Poisson case, otherwise the bounds are much larger. In a slightly different 
direction, Cai and Willmot 0, Theorem 3.2] showed that if {Q(i)} is decreasing then the cumulative 
distribution function of the compound Poisson distribution CPo(A, Q), evaluated at the integers, is 
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log-concave. Finally, Keilson and Sumita 371 Theorem 4.9] proved that, if Q is log-concave, then the 
ratio, 

C Q U x (n) 
C Q n A (n + l)' 

is decreasing in A for any fixed n. 

In the present context, we first show that a compound Bernoulli sum is log-concave if the compounding 
distribution Q is log-concave and the Bernoulli parameters are sufficiently large. 

Lemma 5.1. Suppose Q is a log-concave distribution on N, and all the elements pi of the parameter 
vector p = (pi,p2, ■ ■ ■ ,p n ) satisfy pi > Y+Q^rp/Q^y- Then the compound Bernoulli sum distribution 
Cqbp is log-concave. 

Proof. Observe that, given that Q is log-concave, the compound Bernoulli distribution CBern (p, Q) is 
log-concave if and only if, 

Indeed, let Y have distribution CBern (p, Q). Since Q is log-concave itself, the log-concavity of CBern (p, Q) 
is equivalent to the inequality, Pr(F = l) 2 > Pr(Y = 2) Pr(y = 0), which states that, (pQ(l)) 2 > 
(1 — p)pQ(2), and this is exactly the assumption (|30]) . 

The assertion of the lemma now follows since the sum of independent log-concave random variables is 
log-concave; see, e.g., [3]. □ 

Next we examine conditions under which a compound Poisson measure is log-concave, starting with a 
simple necessary condition. 

Lemma 5.2. A necessary condition for CPo(A, Q) to be log-concave is that, 

A> 2g(2) r<m 

A -w (31) 

Proof. For any distribution P, considering the difference, CqP(1) 2 — CqP(0)CqP(2), shows that a 
necessary condition for CqP to be log-concave is that, 

(P(l) 2 - P(0)P(2))/P(0)P(1) > Q(2)/Q(l) 2 . (32) 

Now take P to be the Po(A) distribution. □ 

Similarly, for P = b pi a necessary condition for the compound Bernoulli sum Cq6 p to be log-concave is 
that, 

Pi , (sr* Pi Uv- ft 2( 9( 2 ) 



^1-Pi + 1^(1^) 2 J t?i-ft 



g(i) 2 ' 



which, since the left-hand-side is greater than J2iPi/(^ ~ Pi) — TliPi> hold as long as J^iPi — 
2Q(2)/Q(1) 2 . 
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Note that, unlike for the Poisson distribution, it is not the case that every compound Poisson distribution 
CPo(A, Q) is log-concave. 

Next we show that for some particular choices of Q and general compound distributions CqP, the above 
necessary condition is sufficient for log-concavity. 

Theorem 5.3. Let Q be a geometric distribution on N. Then CqP is log-concave for any distribution 
P which is log-concave and satisfies the condition (|32p . 

Proof. If Q is geometric with mean 1/a, then, Q* y (x) = a y (l — a) x ~ y (yZ]) > which implies that, 

X 

c Q p( x ) = '£p(y)ay(i-ay-y 

Condition (j3"2"j) ensures that CqP(1) 2 — CqP(0)CqP(2) > 0, so, taking z = y — 1, we need only prove 
that the sequence, 

X 

C{x) := C Q P(x + 1)/(1 - ocf = £ P(z + 1) 

z=0 

is log-concave. However, this follows immediately from [34], Theorem 7.3], which proves that if {a^} is 
a log-concave sequence, then so is {bi}, defined by bi = J2]=o \p a r ^ 

Theorem 5.4. Let Q be a distribution supported on the set {1,2}. Then the distribution CqP is 
log-concave for any ultra-log- concave distribution P with support on {0,1,..., N} (where N may be 
infinite), which satisfies 

(x + l)P(x + l)/P(x)>2Q(2)/Q(l) 2 (33) 

for allx = 0,l,...,N. 

In particular, if Q is supported on {1,2}, the compound Poisson distribution CPo(A, Q) is log-concave 
for allX>^$. 

Note that the condition fl33|) is equivalent to requiring that NP(N)/P(N - 1) > 2Q(2)/Q(1) 2 if N is 
finite, or that lim^o^x + l)P(x + l)/P(x) > 2Q(2)/Q(1) 2 if N is infinite. 

The proof of Theorem 15.41 is based in part on some of the ideas in Johnson and Goldschmidt 29] , and 
also in Wang and Yeh [581 ] . where transformations that preserve log-concavity are studied. Since the 
proof is slightly involved and the compound Poisson part of the theorem is superseded by Theorem [53] 
below, we give it in the appendix. 

Lemma l5~2l and Theorems l5.3l and l5.4[ supplemented by some calculations of the quantities CqIIa(2;) 2 — 
CqI1\(x — 1)CqIIa(ie + 1) for small x, suggest that compound Poisson measure CPo{\,Q) should be 
log-concave, as long as Q is log-concave and AQ(1) 2 > 2Q(2). Indeed, the following slightly more 
general result holds; see Section 7 for some remarks on its history. As per Definition 13.31 we use Q# 
to denote the size-biased version of Q. Observe that log-concavity of Q# is a weaker requirement than 
log-concavity of Q. 
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Theorem 5.5. If is log-concave and XQ(1) 2 > 2Q(2) with Q(l) > 0, then the compound Poisson 
measure CPo(A, Q) is log-concave. 

Proof. It is well-known that compound Poisson probability mass functions obey a recursion formula: 

k 

kC Q H x (k) = \J2jQU)C Q IL x (k - j) for all k G N. (34) 

(This formula, which is easy to prove for instance using characteristic functions, has been repeatedly 
rediscovered; the earliest reference we could find was to a 1958 note of Katti and Gurland mentioned 
by N. de Pril 16 1, but later references are Katti [IBj], Adelson [l| and Panjer (4||; in actuarial circles, 
the above is known as the Panjer recursion formula.) For notational convenience, we write [Iq for the 
mean of Q, rj = X(j + l)Q(j + 1) = A^qQ#(j), and Pj — Cgll\(j) for j £ Z + . Then (jM]) reads, 

fe 

(k + l)p k +l = r lPk-] 
3=0 

for all k G Z+. 

Theorem [53] is just a restatement using (|3"4l of [2^, Theorem 1]. For completeness, we sketch the proof 
of Hansen [22.], which proceeds by induction. Note that one only needs to prove the following statement: 
If Q# is strictly log-concave and XQ(1) 2 > 2Q(2), then the compound Poisson measure CPo(A, Q) is 
strictly log-concave. The general case follows by taking limits. 

By assumption, \Q(1) 2 > 2Q(2), which can be rewritten as Tq > r*i, and hence, 

2 

P 2 i-PoP2 = y^-rOX). 

This initializes the induction procedure by showing that the subsequence (po,Pi,P2) is strictly log- 
concave. Hansen 22] developed the following identity, which can be verified using the recursion (|34l) : 
Setting p_i = 0, 

m(m + 2)[p 2 m+l - p m p m+2 ] = p m +i(rop m - p m +\) 

Av^ (35) 
+ 2_^(p m -iPm-k-i - Pm-kPm-l-i)(rk+ifi - r l+ ir k ). 

1=0 k=0 

Observe that each term in the double sum is positive as a consequence of the induction hypothesis 
(namely, that the subsequence (po,pi, . . . ,p m +i) is strictly log-concave), and the strict log-concavity of 
r. To see that the first term is also positive, note that the induction hypothesis implies that pk+i/pk is 
decreasing for k < m + 1; hence, 

Pi ^ Pm+l 

r = — > . 

Po Pm 

Thus it is shown that p^+i > PmPm+2, which proves the theorem. □ 
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We note that Hansen's remarkable identity (1551) is reminiscent of (although more complicated than) an 
identity that can be used to prove the well-known fact that the convolution of two log-concave sequences 
is log-concave. Indeed, as shown for instance in Liggett 0, if c= a-kb, then, 

c 2 k - Ck-iCk+i = ^(a,aj_i - ai-iaj)(bk-ibk-j+i - bk-i+ibk-j). 

i<j 

Observe that ([34]) can be interpreted as saying that the size-biased version of CqI1\ is the convolution 
of the sequence r with the sequence p. 

6. Applications to Combinatorics 

There are numerous examples of ultra-log-concave sequences in discrete mathematics, and also many 
examples of interesting sequences where ultra-log-concavity is conjectured. The above maximum entropy 
results for ultra-log-concave probability distributions on Z + yield bounds on the "spread" of such ultra- 
log-concave sequences, as measured by entropy. Two particular examples are considered below. 

6.1. Counting independent sets in a claw-free graph 

Recall that for a graph G — (V, E) , an independent set is a subset of the vertex set V such that no 
two elements of the subset are connected by an edge in E. The collection of independent sets of G is 
denoted 1(G). 

Consider a graph G on a randomly weighted ground set, i.e., associate with each i 6 V the random 
weight Xi drawn from a probability distribution Q on N, and suppose the weights {Xi : i £ V} are 
independent. Then for any independent set I G 1(G), its weight is given by the sum of the weights of 
its elements, 

= ]Tx;, 

iel i=l 

where = denotes equality in distribution, and X[ are i.i.d. random variables drawn from Q. For the 
weight of a random independent set I (picked uniformly at random from 1(G)), one similarly has, 

III 

iel »=i 

and the latter, by definition, has distribution CqP, where P is the probability distribution on Z + 
induced by |I|. 

The following result of Hamidoune [2lJ (see also Chudnovsky and Seymour [l3| for a generalization 
and different proof) connects this discussion with ultra-log-concavity. Recall that a claw-free graph is 
a graph that does not contain the complete bipartite graph ^ as an induced subgraph. 

Theorem 6.1 (Hamidoune (2l|). For a claw-free finite graph G, the sequence {Ik}, where Ik is the 
number of independent sets of size k in G, is ultra-log-concave. 
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Clearly, Theorem 16. II may be restated as follows: For a random independent set I, 

P(k) :=Pr{|I| =k}cKl k , 

is an ultra-log-concave probability distribution. This yields the following corollary. 

Corollary 6.2. Suppose the graph G on the ground set V is claw-free. Let I be a random independent 
set, and let the average cardinality of I be A. Suppose the elements of the ground set are given i.i.d. 
weights drawn from a probability distribution Q on N, where Q is log-concave with Q(l) > and 
\Q(1) 2 > 2Q(2). IfW is the random weight assigned to I, then, 

H(W) < H(C Q U X ). 



If Q is the unit mass at 1, then W = |I|, and Corollarv l6.2| gives a bound on the entropy of the cardinality 
of a random independent set in a claw-free graph. That is, 

H(\l\) < H(Tl\), 

where A = E\l\. Observe that this bound is independent of n and depends only on the average size of a 
random independent set, which suggests that it could be of utility in studying sequences associated with 
graphs on large ground sets. And, although the entropy of a Poisson (or compound Poisson) measure 
cannot easily expressed in closed form, there are various simple bounds [Til Theorem 8.6.5] such as, 



H(U X ) < - log 



27re A 



1 

12 



(36) 



as well as good approximations for large A; sec, e.g., [4(| 26, 17 1. One way to use this bound is via the 
following crude relaxation: Bound the average size A of a random independent set by the independence 
number a(G) of G, which is defined as the size of a largest independent set of G. Then, 



H(\l\)< ^ log 



2ne[a(G) + - 



(37) 



which can clearly be much tighter than the trivial bound, < loga(G), using the uniform distri- 

bution, when a(G) > 16. 



6.2. Mason's conjecture 

Recall that a matroid M on a finite ground set [n] is a collection of subsets of [n] , called "independent 
sets"Q, satisfying the following: (i) The empty set is independent, (ii) Every subset of an independent 
set is independent, (iii) If A and B are two independent sets and A has more elements than B, then 
there exists an element in A which is not in B and when added to B still gives an independent set. 

Consider a matroid M on a randomly weighted ground set, i.e., associate with each i £ [n] the random 
weight Xi drawn from a probability distribution Q on N, and suppose the weights {Xi : i G [n]} are 



4 Note that although graphs have associated cycle matroids, there is no connection between independent sets of matroids 
and independent sets of graphs; indeed, the latter are often called "stable sets" in the matroid literature to distinguish 
the two. 
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independent. As before, for any independent set / € M, its weight is given by the sum of the weights 
of its elements, and the weight of a random independent set I (picked uniformly at random from M), 
is, 

iSl i=l 

where the X[ are i.i.d. random variables drawn from Q. Then w(T) has distribution CqP, where P is 
the probability distribution on Z + induced by |I|. 

Conjecture 6.3 (Mason [46|). The sequence {Ik}, where Ik is the number of independent sets of size 
k in a matroid on a finite ground set, is ultra-log-concave. 



Strictly speaking, Mason's original conjecture asserts ultra-log-concavity of some finite order (not defined 
in this paper) whereas this paper is only concerned with ultra-log-concavity of order infinity; however 
the slightly weaker form of the conjecture stated here is still difficult and open. The only special cases 
in which Conjecture 16.31 is known to be true is for matroids whose rank (i.e., cardinality of the largest 
independent set) is 6 or smaller (as proved by Zhao [60]), and for matroids on a ground set of If or 



smaller (as proved by Kahn and Neiman |33J). Very recently, Lenz [42J proved that the sequence {Ik} 
is strictly log-concave, which is clearly a weak form of Conjecture 16.31 

Conjecture 16.31 equivalently says that, for a random independent set I, the distribution, P(k) = 
Pr{|I| = k} oc Ik, is ultra-log-concave. This yields the following corollary. 

Corollary 6.4. Suppose the matroid M on the ground set [n] satisfies Mason's conjecture. Let I be a 
random independent set of M, and let the average cardinality of I be X. Suppose the elements of the 
ground set are given i.i.d. weights drawn from a probability distribution Q on N, where Q is log-concave 
and satisfies Q(l) > and AQ(1) 2 > 2Q(2). IfW is the random weight assigned to 1, then, 

H(W) < H(CqTL\). 



Of course, if Q is the unit mass at f , Corollary 16.41 gives (modulo Mason's conjecture) a bound on the 
entropy of the cardinality of a random independent set in a matroid. That is, 

H(\I\) < H(U X ), 

where A = E\I\. As in the case of graphs, this bound is independent of n and can be estimated in terms 
of the average size of a random independent set (and hence, more loosely, in terms of the matroid rank) 
using the Poisson entropy bound (|36[) . 



7. Extensions and Conclusions 



The main results in this paper describe the solution of a discrete entropy maximization problem, under 
both shape constraints involving log-concavity and constraints on the mean. Different entropy problems 
involving log-concavity of continuous densities have also been studied by Cover and Zhang |f5[ and by 
Bobkov and Madiman Q, using different methods and motivated by different questions than those in 
this work. 
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The primary motivation for this work was the development of an information-theoretic approach to 
discrete limit laws, and specifically those corresponding to compound Poisson limits. Above we have 
shown that, under appropriate conditions, compound Poisson distributions have maximum entropy 
within a natural class. This is analogous to the maximum entropy property of the Gaussian and Poisson 
measures, and their corresponding roles in Gaussian and Poisson approximation, respectively. Moreover, 
the techniques introduced here - especially the introduction and analysis of a new score function that 
naturally connects with the compound Poisson family - turn out to play a central role in the development 
of an information-theoretic picture of compound Poisson limit theorems and approximation bounds Q ■ 

After a preliminary version of this paper was made publicly available [30j j , Y. Yu (59j provided different 
proofs of our Theorems 12.41 and 12.51 under less restrictive conditions, and using a completely different 
mathematical approach. Also, in the first version of (30j . motivated in part by the results of Lemma l5.2l 
and Theorems l5.3l and l5.4[ we conjectured that the compound Poisson measure CPo(A, Q) is log-concave, 
if Q is log-concave and AQ(1) 2 > 2Q(2). Y. Yu [59[ subsequently established the truth of the conjecture 
by pointing out that it could be proved by an application of the results of Hansen in 22]. Theorem [53] 
in Section 5 is a slightly more general version of that earlier conjecture. Note that in order to prove 
the conjecture it is not necessary to reduce the problem to the strictly log-concave case (as done in the 
proof of Theorem 15. 5p , because the log-concavity of Q implies a bit more than log-concavity for Q# . 
Indeed, the following variant of Theorem 15.51 is easily proved: If Q is log-concave with Q(l) > and 
AQ(1) 2 > 2Q(2), then the compound Poisson measure CPo(A,Q) is strictly log-concave. 

In closing we mention another possible direction in which the present results may be extended. Suppose 
that the compounding distribution Q in the setup described in Section 2 is supported on R and has a 
density with respect to Lebesgue measure. The definition of compound distributions CqP (including 
compound Poissons) continues to make sense for probability distributions P on the nonnegative integers, 
but these now clearly are of mixed type, with a continuous component and an atom at 0. Furthermore, 
limit laws for sums converging to such mixed-type compound Poisson distributions hold exactly as in the 
discrete case. It is natural and interesting to ask for such 'continuous' analogs of the present maximum 
entropy results, particularly as neither their form nor method of proof are obvious in this case. 
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Appendix A. Appendix 



Proof. Writing R(y) = y\P(y), we know that C Q P(x) = YT v =o R (v) (Q* v (x)/y\) . Hence, the log- 
concavity of CqP(x) is equivalent to showing that, 

V V R{y)R{z) (r) (Q^pO^M _ Q* v {x + l)Q* z {x-l)\ 

for all x > 2, since the case of x — 1 was dealt with previously by equation (1321) . In particular, for (i), 
taking P = Po(A), it suffices to show that for all r and x, the function, 



«(*):= £ 

y-\-z— r 



r\ Q* y {k)Q* z (2x-k) 



vj Q* r (2x) 



is unimodal as a function of k (since g r ,x{k) is symmetric about x). 

In the general case (ii), writing Q(2) = p = 1 - Q(l), we have, Q* y (x) = ( x v _ y )p x ~ y (l —p) 2y ~ x , so that, 

'r\ Q* y (k)Q* z (2x -k) _ (2x - r\ (2r - 2x\ 

(A. 2) 



y) Q* r {2x) \k-y J \2y - k 

for any p. Now, following [29l . Lemma 2.4] and (58l . Lemma 2.1], we use summation by parts to show 
that the inner sum in (|A.ip is positive for each r (except for r — x when x is odd), by case-splitting 
according to the parity of r. 

(a) For r = 2t, we rewrite the inner sum of equation (|A.1|) as, 
t 

^2(R(t + s)R(t - s) - R(t + s + l)R(t -s- 1)) x 



s=0 

/ t+s 



m-^ I / 2x - r\ f 2r - 2x\ f 2x - r \ f 2r - 2x 
\\x-y)\2y-x) \x + 1 - y) \2y - x - 1 



\y=t-s 



where the first term in the above product is positive by the ultra-log-concavity of P (and hence log- 
concavity of R), and the second term is positive by Lemma | Appendix A.l below. 

(b) Similarly, for x ^ r = 2t + 1, we rewrite the inner sum of equation (|A.1[) as, 
t 

^2(R(t + s + l)R(t-s)-R(t + s + 2)R(t-s-l))x 



s=0 

/t+l+s 

1 E 



\v=t 



2x - r\ (2r - 2x\ ( 2x - r \ ( 2r - 2x 
x-y )\2y-x ) \x + l~yj \2y -x-1 



where the first term in the product is positive by the ultra-log-concavity of P (and hence log-concavity 
of R) and the second term is positive by Lemma | Appendix A.l| below. 
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(c) Finally, in the case of x = r = 2t + 1, substituting k — x and k = x + 1 in (|A.2[) . combining the 
resulting expression with (|A.1[) . and noting that ( 2r ~ 2a: ) is 1 if and only if u = (and is zero, otherwise), 
we see that the inner sum becomes, —R(t + l)i?(i)( 2t J 1 " 1 ), and the summands in (|A.1[) reduce to, 

p x R(t)R(t + 1) 
(t + l)W. • 

However, the next term in the outer sum of equation (|A.1|) . r = x + 1, gives 



> 



2(2i)! 

p*' 1 ^ -p) 2 
2(2*)! 



R(t + 1) 2 \ 2 
R(t + 1) 2 



2t\ ( 2t 
t)~\t+l 
2t\ ( It 
t J [t + l 



R(t)R(t + 2) 



2t 



p x - 1 {\-pY 

2{t + \)\t\ 



R(t + 1) 



Hence, the sum of the first two terms is positive (and hence the whole sum is positive) if R(t + 1)(1 — 
p) 2 /{2p) > R(t). 

If P is Poisson(A), this simply reduces to equation (|31l) . otherwise we use the fact that R(x + 1)/R(x) 
is decreasing. □ 

Lemma Appendix A.l. (a) If r — 2t, for any < s < t, the sum, 

t+s 



E 

V=t—s 



2x-r\(2r~ 2x 
x-y )\2y-x 



2x-r \f 2r~2x 
x + 1 — yj \2y — x — 1 



> 0. 



(b) If x r = 2t + 1, for any < s < t, the sum, 

t+l+s 



E 

y=t—t 



2x — r 
x-y 



2r - 2x 
2y-x 



2x ~ r 
x + 1 — y 



2r - 2x 
2y-x-l 



> 0. 



Proof. The proof is in two stages; first we show that the sum is positive for s = t, then we show that 
there exists some S such that, as s increases, the increments are positive for s < S and negative for 
s > S. The result then follows, as in [29[ or (58j . 



For both (a) and (b), note that for s — t, equation (|A.2|) implies that the sum is the difference between 
the coefficients of T x and T x+1 in f r , x {T) = (1 + T 2 ) 2x - r {l + T) 2r - 2x . Since f r ,x{T) has degree 2x 
and has coefficients which are symmetric about T x , it is enough to show that the coefficients form a 
unimodal sequence. Now, (1 + j , ^2x-r^^ _|_ coefficients which do form a unimodal sequence. 

Statement Si of Keilson and Gerber [361 ] states that any binomial distribution is strongly unimodal, 
which means that it preserves unimodality on convolution. This means that (1 + x 2 ) 2x ~ r (l + T) 2r ~ 2x 
is unimodal if r — x > 1, and we need only check the case r = x, when f r ,x{T) = (1 + T 2 ) r . Note that 
if r = 2t is even, the difference between the coefficients of T x and T x+1 is ( 2 t *), which is positive. 

In part (a), the increments are equal to ( 2 ^7+* ) i2t-~2s-x) multiplied by the expression, 



(x-t- s)(2t -2s- . 



(x - t + s)(2t + 2s - x) 



(x + 1 - t + s)(2t + 2s - x + 1) (x + l-t- s)(2t - 2s - x + 1) 



2G 



which is positive for s small and negative for s large, since placing the term in brackets over a common 
denominator, the numerator is of the form (a — bs 2 ). 

Similarly, in part (b), the increments equal (^Ijl" 1 ) (^^-x) times the expression, 



2 - 



x-t+s ) \2t-2s-xJ 

(x-t-s- l)(2t -2s-x) (x-t + s)(2t + 2 + 2s - x) 



(x + 1 - t + s)(2t + 2s - x + 3) (x-t- s)(2t + 1 - 2s - x) ' 
which is again positive for s small and negative for s large. □ 
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